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^ We report on a lattice QCD calculation with two dynamical flavors of the isovector vector corre- 

lator in the high- temperature phase. We analyze the correlator in terms of the associated spectral 
\^ function by performing a fit for the difference of the thermal and vacuum spectral functions, using 

also an exact sum rule that constrains this difference. Additonally we carry out a direct fit for the 
thermal spectral function, and obtain good agreement between the two analyses for frequencies 
below the two-pion threshold. Under the assumption that the spectral function is smooth in that 
region, we give an estimate of the electrical conductivity. 
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1. Introduction 



The properties of strongly interacting matter under extreme conditions are the subject of in- 
tensive experimental and theoretical investigation. A comprehensive picture of a state of matter 
requires not only the knowledge of equilibrium properties such as the equation of state and static 
susceptibilities, but also an understanding of its transport properties. 

Here we report on recent progress in determining the Euclidean isovector vector correlation 
function using dynamical light quarks in the high temperature phase of QCD. This enables the first 
analysis of its underlying spectral function with unquenched quarks. We are able to determine the 
gross features of the thermal spectral function by analyzing directly the difference of the thermal 
and vacuum correlators. This difference can be further constrained using a recently derived sum 
rule [1]. In addition we employ an approach followed recently in a quenched study to determine 
the spectral function directly from the thermal correlator based on an appropriate Ansatz [2-4]. 

Note, this proceedings article represents a shortened version of [5] and we refer the reader to 
this publication for a more in depth discussion. 



2. Basic definitions and expectations 

Our primary observables are the Euclidean vector current correlators and their spectral repre- 
sentation: 

. J.'. MO)', ^ /; g a, 

with/^(x) = -^[ii{x)Ynu{x) —d{x)Y^d{x)) the isospin current. For a given function p(a), T), the 
reconstructed correlator is defined as 

Jo Ik smh(G)p/2) 

It can be interpreted as the Euclidean correlator that would be realized at temperature T if the 
spectral function was unchanged between temperature T and T' . For T' = it can be directly 
obtained from the zero-temperature Euclidean correlator via [6] 

G'^\xJ) = G'''''{tJ;Q)= £ G(|T + mj3|,r = 0). (2.3) 



In the thermodynamic limit, the subtracted vector spectral function obeys a sum rule (see [1] sec. 
3.2), 

d(0 

/ — Ap(w,r) = o, Ap(w,r) = p,7(w,r)-p,7(w,o). (2.4) 

J-oo ft) 



The electrical conductivity of QCD, connected to the isospin diffusion constant D, is given by 
a Kubo formula in terms of the low-frequency behavior of the spectral function, (where Cem = 

i:/=Me} = 5/9) 

a = C„„D;t. = %lim^^^^^. (2.5) 
6 ffl^o CO 
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64^ X 128 


64^ X 16 


Ref. 




64^ X 128 


64^ X 16 




5.50 




T [MeV] 




253(4) 


K 


0.13671 








0.871(1) 




1.7515 






4.42(31) 




Zv 


0.768(5) 


[7] 


mi/T 


3.33(5) 




^[fm] 


0.0486(4)(5) 


[8] 




1.244(5) 




m;t[MeV] 


270 


[8] 


a/T 


5.98(11) 





Table 1: The left part of the table shows the common quantities characterizing the zero-temperature and 
finite-temperature ensembles. In the right part, the fit parameters for the vacuum correlator in units of 
T = 253MeV and the value of the (isospin) quark number susceptibility Xs/T^ cire given. For more details 
on the generation of the = 128 ensemble, see [8]. The number of configurations generated with Nx = 16 
is 317. 

3. Lattice QCD data 

All our numerical results were computed on dynamical gauge configurations with two light, 
mass-degenerate ^(a)-improved Wilson quark flavors. We calculated correlation functions us- 
ing the same discretization and masses as in the sea sector in two different ensembles. The first 
corresponds to virtually zero-temperature on a 64^ x 128 lattice (labeled 07 in [8]) with a lattice 
spacing of a = 0.0486(4) (5)fm [8] and a pion mass of nij^ = 270MeV, so that mj^L = 4.2. Sec- 
ondly we generated an ensemble on a lattice of size 64^ x 16 with all bare parameters identical 
to the zero-temperature ensemble. In this way it is straightforward to compare the coiTclation 
functions respectively in the confined and deconfined phases of QCD. Choosing A^^ = 16 yields a 
temperature of T ~ 250MeV. Based on preliminary results on the pseudo-critical temperature 
of the crossover from the hadronic to the high- temperature phase [9], the temperature can also be 
expressed asT /T^^ 1.2. 

The vacuum correlator serves as a reference in this work. To fix the parameters of the lightest 
vector state in the finite volume of the simulation, we fitted the vacuum correlation function to an 
Ansatz of the form^ 

G;,(t,0) =Aie-"''' + -^Ko exp(-nT) (O? /t + 2^/^^- + 2/t^) . (3.1) 

In addition we estimate the thermal (isovector) quark number susceptibility Xs/T'^ from the 
time-time component of the vector coiTelation function. The parameters used in our lattice setup, 
the 'p -meson' parameters and the value of the static susceptibility are summarized in Tab. 1. 

3.1 Thermal and vacuum correlators 

In Fig. l(left) we show the correlator Gii{T,T) computed at T ~ 250MeV together with the 
corresponding free 'continuum' and free 'lattice discretized' correlation functions. In addition 
we show the reconstructed correlator G^f'^{z) as obtained from Eq. 2.3. The reconstructed cor- 
relator lies somewhat lower than the thermal correlator. The insert in Fig. 1 displays the ratio 
Gii{z) / G^f^ (z) in order to make their relative r dependence visible. For small Euclidean times 

'For details please see [5] 



3 



Two-flavour lattice QCD correlation functions in the deconflnement transition region 



Anthony Francis 



1e+04 5 ^ ^ ^ ^ ^ ^ ^ ^ 1 




12345678 12345678 



Figure 1: Left: Thermal Gu{z) /T^ and reconstructed G[J^'^(t) /T^ vector correlators over Euclidean time 
separation T compared to the free ( continuum and lattice) cases. The reconstructed correlator was computed 
by applying Eq. 2.3 to the data obtained from a lattice sized Na ~ 64 and = 128 . The insertion shows the 
ratio G,-,(t) /GJ" (t). Right: Difference AG(t, T)/T^ of the thermal vector correlator and the corresponding 
reconstructed correlator as a function of Euclidean time T. The insertion shows AG(t, T)/T^ in the region 
x/a > 3. 



T < jS /4 this ratio is unity, above it increases monotonically until it levels off around the midpoint 
at about 10% above unity. A thermal modification of the spectral function has thus taken place (re- 
call that the spectral function underlying GJJ^ (t) contains the bound states of the confined theory). 
In Fig. 1 (right) we show the difference 

AG{T,T) = Gii{i:,T)-Gij {t:,T)= — Ap(w,r — .,,n/^x — (3-2) 

of the thermal and the reconstructed correlators. Given that p;,(G),r) and p;,(ft),r = 0) have the 
same ~ co^ behavior, this means we are subtracting non-perturbatively the ultraviolet tail of the 
spectral function. Using this difference we are therefore able to probe the change in the vector 
spectral function from the confined to the deconfined phase for frequencies (O < 0{T). 

For small times, the difference (3.2) turns out to be negative, while it is positive for T > j3 /4. 
Note the en^ors decrease with increasing Euclidean time throughout the available range. We show 
a more detailed view of the region r/a > 3 in the insert of Fig. 1 (right). Here the difference still 
exhibits a mild increase and levels off near the midpoint. The value it reaches at the midpoint is 
AG(t = j8 /2, r) /T^ = 0.291 (55). 



4. Analysis of lattice correlators in terms of spectral functions 

4.1 Fit to the thermal part of the vector correlator 

We proceed to investigate the behavior of the thermal part of the spectral function Ap by 
fitting the difference of the thermal and the reconstructed correlator, see Eq. (3.2). As described 
in the previous subsection, the fact that the data (displayed in Fig. 1 (right)) is positive at long 
distances and negative at short distances suggests that the thermal spectral weight exceeds the 
vacuum spectral weight at low frequencies and falls short of it at higher frequencies. 
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Figure 2: Left: Fits to AG{t,T)/T^ = [G,i(T) - G;7'(T)]/rl The blue and red results differ by the form 
of transport peak in the Ansatz- The error bands are computed from the covariance matrix of the flt. Right: 
The resulting spectral functions for both Ansdtze. 

We thus parametrize Ap using the following Ansatz for ft) > 0: 

Ap(w,r) = pT{co,T)-pB{co,T)+^pF{(0,T), (4.1) 
IcBgB tanh(co/r)3 

Pb{(0,T) = —— , 2 ' (4.2) 

4{(0-mBy + gi 

/ 4cw , ^, 4cTtanh((o/T) , ^ 

(a)/g)2 + l (co/g)^ + l 



ApFi(0,T) = Pf{(O,T)-Pf{(O,0), pf(w,r) = ^JCCo2tanh(^) 



(4.4) 



The bound state (B) and the transport peak (T) are represented by Breit-Wigner forms. Even such 
a simple Ansatz requires three parameters {cB,gB,mB) to determine the bound state peak, two 
parameters {c,g) for the transport peak and one {k) for the 'perturbative' contribution (F). We will 
therefore fix some of them using the vacuum correlator. In the following we set equal to mi, 
given in Tab. 1 , which we obtained from the exponential fit to the vacuum correlator. Note that the 
ai^ea under the bound state peak / dw Pb/(0 does not depend on the width gs in the limit where it is 
small. We therefore perform fits for three fixed values of this parameter, and check the sensitivity 
of the result. We choose the values gg/r = 0.1,0.5 and 1.0. 

The tail ~ T/(0 of the Ansatz pr i violates the OPE prediction that Ap ~ {T/co)^ at large 
frequencies. It has been ai^gued in [10] that this might lead to an overestimate of the transport 
contribution. To avoid this problem we introduce the Ansatz 2, where CO — >■ rtanh(ft)/r). This 
Ansatz possesses the correct asymptotic behavior, as well as the expected linear behavior in (O at 
small frequencies. Finally, to complete the parametrization of Ap{(o), we include a weak-coupling 
term describing the subtraction of the large frequency parts of the thermal and vacuum spectral 
functions. This contribution Pf(w, K") — vanishes exponentially as the frequency increases. 

In the next step we fit the combined Ansatze of Ap{(0,CB,gB,mB,c,g, k) to the data, while at 
the same time satisfying the sum rule of Eq. 2.4 to an accuracy of 10^^. We limit ourselves to 
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Figure 3: Left panel: fit results for G,-, (t) /G^l^" {t) using both Ansdtze for the transport peak. The resulting 
ratio of thermal moments /;(2.0)/7;(2 0)/><.c displayed on the right side of the plot. Right panel: The 
corresponding spectral functions normalized by COT. 

fitting tiie region 5 < r/a <S only, in order to minimize the influence of cut-off effects. With mb 
determined by the vacuum correlator, we set gs successively to the three different values mentioned 
above and fixed K around unity, and fitted the parameters c, g and cb- The errors and error bands 
shown in the following have been computed using the covariance matrix of the corresponding fit 
for fixed values of and K. 

The resulting correlators and spectral functions are displayed in Fig. 2 for gB/T = 0.50 and 
K" = 1.10. In the left panel of Fig. 2 the data AG{x,T) is compared to the fits using Pt,i{o)) and 
Pt;i{(o) as transport contribution. We achieve a quasi -perfect description of the data for r/a > 4. 
The right panel shows that both Ansatze exhibit a substantial spectral weight around the origin and 
a negative contribution from the region of the p mass. 

4.2 Weak-coupling inspired fit to the thermal vector 

In contrast to the previous section, here we study directly the thermal vector correlator and its 
ratio of thermal moments /?(^'°', computed along the lines of [2]. We perform a fit inspired by the 
weak coupling form of the thermal spectral function, 

p{(0,T) = Pt{(0,T)+Pf{(0,T) , (4.5) 

where the form of the two contributions is defined in Eq. (4.3) and (4.4). At a given temperature 
this Ansatz is characterized by three parameters {c,g, k). We fit the full Ansatz p{(0,c,g, k) to 
the thermal correlator G,;(t), while at the same time demanding that be reproduced. In this 
analysis the three parameters c,g and K are fitted, and the fit range is5<T/a<8as before. 

The resulting correlators and spectral functions are shown in Fig. 3. The ratio Gii{r)/G^i^^^{T) 
in the left panel of Fig. 3 is well described by both versions j and Pt,2 of the transport contri- 
bution for x/a > 5, while also the ratio of thermal moments (given on the far right of the plot) is 
reproduced. For x/a <5 our Ansatz fails to reproduce these points, which we suspect is partly due 
to cutoff effects. 
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Figure 4: Comparison of the spectral functions obtained from analyzing (a) AG{t,T) and (b) Gii{T,T) 
using in both cases P7-.2(fi)) ~ tanh(a)/r) in the low frequency region. All curves have been multiplied by a 
factor 1/6 and divided by COT, entailing that the intercept at CO ~0 yields an estimate of O /CgmT- 

On the right hand side of Fig. 3 we show the resulting spectral functions divided by (oT . 
Clearly both Ansatze give very similar results that lie within errors of each other. The thermal 
correlator is even less sensitive to the asymptotic behavior of the transport contribution in the 
Ansatz than in the difference of correlators studied in Sec. 4.1. 

5. Discussion 

We now compare the results of the fits to AG(t, T) and G,7(t, T). Since the vacuum spectral 
function vanishes below Inij^ k, 540MeV (in infinite volume), pij{(0,T) and Ap(w,r) should be 
equal for co < Irtin ^ 2. IT. We thus plot the spectral functions obtained from the two fits in 
this frequency region, see Fig. 4. Here we restrict ourselves to showing these results based on 
Pt;2{(0, T), as their theoretical foundation is more sound than those with pr i (ft), T). All curves are 
multiplied by a factor 1/6, which means that the intercept at ftj = yields an estimate of a/Ce,„r 
with o the electrical conductivity of the quai^k gluon plasma. 

The results obtained by fitting Gii{z,T) agree very well with the central values obtained by 
fitting AG{z,T), whereby the fit to AG{z,T) using the transport Ansatz pr2(w) yields a slightly 
lower intercept. If we assume the spectral function to be as smooth around the origin as Fig. 4 
suggests, we obtain the following estimate for the electrical conductivity of the quark gluon plasma 
at r ~ 250MeV, 

7^ = 0.40(12), (5.1) 

where Cgm = llf=u,d Qf- Although obtained under a strong assumption, it is interesting to compai^e 
(5.1) to other lattice results obtained under similar assumptions. The following comparison is made 
with quenched results, since to our knowledge there are no previous dynamical QCD studies. 
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Quenched calculations using staggered fermions and different methods for analyzing the spec- 
tral function obtained a/T = ICemV^ 1] and o /T = 0.4(l)Ce,„ [12] in the temperature region 1.5 < 
T /Tc < 3.0. Also, a recent quenched study using Wilson-Clover fermions in the continuum limit 
obtained 0.33Ce„, < <j/T < ICg,,, at T ~ 1.457^ [2] and T ~ l.irt.[4], with similar non-continuum 
results up to temperatures T ~ 2.98rc[3]. Our results using dynamical Wilson-Clover fermions at 
Nr = 16 are thus completely compatible with the recent quenched results. 
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